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We concisely review semileptonic B decays, focusing on recent progress on both theoret- 
ical and experimental sides. 

1. Introduction 

A precise knowledge of semi-leptonic decays of B mesons brings several advantages 
to flavour physics. They can be studied in the context of the heavy flavour effective 
theory. Leptonic and hadronic contributions factorize, and we have a better con- 
trol of the effects of strong interactions, compared to non-leptonic decays. On the 
experimental side, they are not helicity suppressed as leptonic decays. 

Semileptonic decays contribute considerably to the analysis of the unitarity tri- 
angle. The so-called unitarity clock, the circle around the origin in the p— fj plane [l], 
is proportional to the ratio |V^b/T4(,|, which is most cleanly determined from semi- 
leptonic decays. The parameters \V c b\ normalize the whole unitarity triangle. A 
clean determination of \V c b\ and \V u b\ from tree level processes, that are generally 
assumed not significantly affected by new physics at the current level of precision, 
is a valuable input for other new physics sensitive estimates. 

The inclusive and exclusive determinations of \V c b\ rely on different theoretical 
calculations, each with different (independent) uncertainties. The former employs a 
parton-level calculation of the decay rate as a function of the strong coupling con- 
stant and inverse powers of the b-quark mass uib- The latter hinges on a parameteri- 
zation of the decay form factors using heavy quark symmetry and a non-perturbative 
calculation of its normalization. The inclusive and exclusive experimental measure- 
ments make use of different techniques and have, to a large extent, uncorrelated 
statistical and systematic uncertainties. This independence makes the comparison 
of \V c b\ from inclusive and exclusive decays a powerful test of our understanding 
of semileptonic decays. The values of \V u b\ extracted from exclusive and inclusive 
decays differ by about two standard deviation, and a lot of effort, from both exper- 
imental and theoretical side, has been devoted to solve this long-standing puzzle. 
In this short review we discuss semileptonic B decays, focussing on features where 
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recent progress has been reported^ 
2. Exclusive \V c b\ Extraction 

To extract \V c b\ we analyze the tree level driven B — > D^*'lv weak decays. The 
starting point are their differential ratios, needed for taking into account experi- 
mental cuts. They can be parameterized in terms of Q(u>) and according to 
the formulas 

riv r 2 

^-(B-+DlD) = ^ (m B + m D ) 2 m 3 ^ 2 - I) 3 / 2 \V cb \ 2 (G (oj)) 2 
rlV c 2 

— (B -+ D* ID) = ^{m B - m D .) 2 m%.(L, 2 - l) l l 2 \V cb \ 2 X {u)(F(u)f (I) 

neglecting the charged lepton and neutrino masses. G(oS) and ^(cj) are a lin- 
ear combination of the form factors which parameterize the matrix elements 
< D(v D )\V^\B(v B ) > and < D*(vp, o)\A^\B{vb) >, respectively (a refers to the 
D* polarization). is a known phase space factor. The form factors depend 

on uj = vb ■ the product of the heavy quark velocities vb — Pb/^b and 

= PdM / m DM ■ In the B rest frame, u> corresponds to the energy of D^*' nor- 
malized by its mass, that is to w = E D (*) /m D (*). The values of u are constrained 
by kinematics: w > 1, with largest value w ~ 1.5. 

The experiments allow to fit only the product of the form factors multiplied by 
the CKM parameter, that is \ V c b G{u>)\ or \V c b F{u>)\. Due to the kinematic suppres- 
sion factors, (uj 2 — l) 3 / 2 and (ui 2 — l) 1 ^ 2 , data are taken at u) ^ 1. The results are 
then extrapolated to uj = 1, the so-called nonrecoil point, where vb — v D {,), with 
leptons back to back and the mesons at rest. At the nonrecoil point, heavy quark 
symmetries play a useful role in constraining the form factors; with the additional 
input of \G(1)\ and |.F(1)| coming from theory, a value for \V c b\ can be obtained. The 
main theoretical problem is the non perturbative evaluation of the operator matrix 
elements. Lattice and QCD sum rules are the two more common routes toward such 
evaluation. 

Lattice determinations for B — > I v decays are generally more precise than 
in most leptonic and semileptonic decays. Form factors can be described to high 
accuracy by a normalization and a slope. It is possible to connect \V c b\ with the 
semileptonic form factors via ratios (or double ratios) where most uncertainties 
cancel in the heavy-quark symmetry limit. Such ratios tend to reduce most of the 
normalization uncertainty in the lattice currents, as well as highly correlated sta- 
tistical errors in the numerator and denominator. 

The most recent unquenched lattice calculation for B — > DID adopts a 2 + 1 
dynamical gauge configurations generated by the MILC collaboration, and an action 
with the Fermilab interpretation. After correcting by a factor of 1.007 for QED, it 



a D semileptonic decays are a counterpart of the B ones in the charm sector: for a very recent 
review see e.g. |2| 
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gives [3] 

Q{1) = 1.074 ± 0.024 (2) 



The latest more precise data have come in 2008 and 2009 from Babar [4||5]. Using 
the latest HFAG average [6j, which also includes older Aleph, CLEO and Belle 
measurements, the fit results in 

\V cb \Q{l) = (42.64 ± 1.53) x 10~ 3 (3) 

which can be turned into the following estimate: 

\V cb \ = (39.70 ± 1.42«p ± 0.89 t h) x 10" 3 (4) 

With respect to the previous HFAG determination, the theoretical error has been 
slightly reduced, increasing the modest dominance of the experimental uncertainty. 

An alternative lattice approach consists in calculating the form factor normal- 
ization directly at values us > 1, that may allow more precise determinations, 
avoiding the large extrapolation to u) = 1 and reducing model dependence. It is 
currently available only in the quenched approximation, that neglects u, d, and s 
quark loops [7|[8j. This approach, by using 2009 Babar data {5], gives the estimate 

\V ch \ = (41.6 ± 1.8 ± 1.4 ± 0.7 FF ) x 10~ 3 (5) 

a slightly higher value than Q. The errors are respectively statistical, systematic 
and due to the theoretical uncertainty in the form factor G(l). 

Measuring the differential rate for B — > D*lD is easier than for B — > DW since 
the rate is higher and there is no background from mis-reconstructed B — ¥ D*W. 
The most recent HFAG experimental fit f6~| gives 

I Vcfol I J^(l) I = (35.90 ±0.45) x 10" 3 (6) 

while one can obtain [9|[l0] 

T(l) = 0.908 ±0.017 (7) 

by using lattice with the Fermilab action for b— and c— quarks, the asqtad staggered 
action for light valence quarks, and the MILC ensembles for gluons and light quarks. 
It includes the enhancement factor 1.007, due to the electroweak corrections to the 
four-fermion operator mediating the semileptonic decay. This is the latest lattice 
result and updates, one year later, the first (2±l)-flavor lattice calculation of .F(l) [9j 



11 , reducing the total uncertainty on .F(l) from 2.6% to 1.7%. No lattice calculation 
is still available at uj 7^ 1 for B — > D* I v decays. 

By taking Eq. ([6| and Eq. Q, the latest value of |V c b| from exclusive B — > D*lv 
reads [6] 

\V ch \ = (39.54 ± 0.50 cxp ± 0.74th) x 10~ 3 (8) 

where the errors come from experiment and QCD lattice calculation, respectively. 
This number is in excellent agreement with the result Q coming from B — > Dlu 
decay. 
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Non-lattice results generally give one or two sigma lower form factors, and thus 
a larger value of \V c b\ exclusive. The .F(l) form factor has recently been calculated 
using zero recoil sum rules, yielding to 



including full a s and up to 1/ rrc| 



^(1) 
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0.86 ±0.02 



(9) 



Let us remark that in order to compare the 
value in Eq. (|9| with the lattice result Eq. ffin, one has to remove the electroweak 
factor 1.007 from the latter. By using Eq. ([6]) in combination with Eq. ([9]), one 
obtains 



\Vcb\ = (41.6 ±0.6 exp ± 1.9th) x 10" 



(10) 



The precision of non-lattice calculations are affected by deviations from the 
heavy mass limit, which may impact also lattice calculations, as discussed e.g in 
Refe. 13 14 . The deviations are estimated to be parametrically larger in the case 
of B ~ > D with respect to B — s- D* . In order to get a more precise prediction for the 
form factor 5(1), it has been suggested [15] to combine the heavy quark expansion 
with a "BPS" expansion. The latter is an expansion in the limit where the kinetic 
energy /ij is equal to the chromomagnetic moment fi^: in this limit 

5(1) = 1.04 ±0.02 (11) 

With such estimate one finds fl6l 



\V cb \ = (40.7±1.5 cxp ±0.8 th ) x 10 



(12) 



The larger experimental uncertainty with respect to Eq. ( 10 ) reflects the above 
mentioned experimental difficulty in measuring the differential rate for B — » DID. 



3. B— Mesons Decays to Excited D— Meson States 

Let J be the total angular momentum of an heavy meson. In the limit where the 
heavy quark mass is infinity, the spin s of the heavy quark is conserved and decou- 
ples from the total angular momentum of the light degrees of freedom ji = J — s. 
By definition, ji becomes a conserved quantity as well, and its eigenvalue ji a good 
quantum number, despite the fact that the light-degree of freedom consists of a 
superposition of the light quark, sea quarks and gluons. In the non-relativistic con- 
stituent quark model, the open charm system can be classified as n 2S+1 Lj , where n 
is the radial quantum number, L corresponds to the eigenvalue of the relative angu- 
lar momentum between the c-quark and the light degrees of freedom, and S = s®si, 
si being the spin of the light valence quark and s the eigenvalue of the heavy quark 
spin. The experimentally observed charmed mesons are associated with the IS, 25, 
IP, and ID states of the meson wave function. 

The D and D* states correspond to the IS wave function; they form the multiplet 
J p = (0~, 1~), respectively, as given by J = ji ® s = © 1, with s = ji = 1/2. 

When L > 0, states are said orbitally excited. The first orbitally excited bound 
state correspond to the IP wave function of the charmed system (L = 1). Since 
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ji = L®si, it is characterized by ji = 1/2 and jt = 3/2. The former gives the doublet 
with J = j t ® s = 1/2 ® 1/2 = © 1, that is J p = (0+, 1+), denoted as D* and L>J, 
respectively. The latter gives the doublet with J — ji®s — i/2®l/2 = 1 © 2, that 
is J p = (1+.2+), corresponding to Di and D\. The Dq, D[, Di and D\ states are 
generically denoted as D**. 

Other excited discovered states are associated with 2S and ID wave function of 
the charmed system. When n, the radial excitation quantum number, is larger than 
1, as in the 2S state, they are said to be radially excited. 

Parity and angular momentum conservation constrain the decays allowed for 
each state, helping to experimentally identify the D** candidates. The ji = 1/2 
states decay through a S-wave to D(*'tt, and they are both expected to be broad 
(large decay widths), while the ji = 3/2 states decay to the same states through a 
D-wave and are expected to be more narrow states (small widths). 

Recent data coming from e + e~ colliders identify the four candidates to D** 
states with D* (2400), D[ (2430), D x (2420) and £>^(2460). The latest observation of 



B — > D**lvi has been reported by Babar 17 . The rates for the D** narrow states 



are in good agreement with the 2005 measurements by D0 [18]; the ones for the 
broad states are in agreement with DELPHI Jl9|, but do not agree with the D[ limit 
of Belle [20]. Measurements also indicate that the relation expected in the heavy 



mass limit [21-23 



r (b -> d** (j, = f\ w\ > r (b -> d** (j, = w^j (13) 

may be violated; this is known as the "1/2 vs 3/2 puzzle". Moreover, there has been 
a long standing problem with the measured semileptonic branching fractions, that 
can be determined with good precision by integrating the differential decays rates 
for B -> £>W I v decays. The sum of the measured exclusive rates is less than the 
inclusive one: explicitly 

B(B+ -> X c l + v) - B(B + -> Dl + v) - B(B + -> D ( * ] l + v) 

-B{B+ -> Z?(*)7r?+^) = (1.45 ± 0.67)% (14) 
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with a similar relation holding for the corresponding B decays 20 

The broadness of J = 1/2 states may be one reason causing the disagreement 
within experiments and with theory, since it has always been quite difficult to disen- 
tangle very broad resonances from continuum, both on theoretical and experimental 
sides. On these premises, it has been suggested to clarify the comparison between 
theory and experiment analyzing states analogous to Dq and D' ly but narrow, in 
particular studying the decay — > D s jir [25] . Other theoretical suggestions to 
ease or solve the previous problems include taking into account an unexpectedly 
large B— decay rate to the first radially excited D'^ 
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In 2010 Babar has found evidence for two new states 27 , which may be identified 
with the 2S states in the quark model picture. Moreover, new data has recently 
become available on semileptonic B— decays to final states containing a D s ' K 
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system, providing information about the poorly known region of hadro nic masses 
above 2.46 GeV/c , that covers radially excited D— meson states 



28 
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4. Exclusive \V u b\ determination 

The analysis of exclusive charmless semileptonic decays is currently employed to 
determine the CKM parameter \V u b\- Here we focus on the channel B — > irlvi, 
where recent progress has been reported from both experimental and theoretical 

rflv 



31 32 . The 



sides. Other interesting channels are B — > ujIv [30] and B 
B — > irli>i decay is the simplest to interpret, as it is affected by a single form factor 
f+(q 2 ) (in the limit of zero leptonic mass) 

dT(B —} ttWi) G 2 f \ P ^ 3 



2\|2 



W\f+(r)\ 



(15) 



dq 2 24tt 3 

where q is the momentum of the leptonic pair and is the momentum of the pion 
in the B meson rest frame. 

Most studied routes to calculate the form factor are once again lattice and light- 
cone QCD sum rules. Theoretical predictions for the form factor split into two: 
predictions for the form factor normalization /+(0) and for the functional form of 
the q 2 dependence. 

By using current lattice QCD methods, the hadronic amplitudes for B — > nWi 
can be calculated quite accurately because there is only a single stable hadron in 
both the initial and final states. The first results based on unquenched simulations 
have been obtained by the Fermilab/MILC collaboration [33] and the HPQCD col- 



laboration 34 , and they are in substantial agreement. Such high-statistics calcula- 



tions have been performed in the kinematic region where the outgoing light hadron 
carries little energy (q 2 > 16 GeV 2 ). At low q 2 , with light hadrons carrying large 
momentum of order 2 GeV, direct simulations require a very fine lattice which is not 
yet accessible in calculations with dynamical fermions. Such fine lattice would also 
be required to simulate heavy quarks; alternately, one can resort to effective heavy 
quark theory. It is also helpful to rely on extrapolations from larger momentum 
transfer q 2 . In Ref. 33 , the b— quark is simulated by using the so-called Fermilab 



heavy-quark method, while the dependence of the form factor from q is parame- 



terized according to the z-expansion 35 -37 . In Ref. [34] , the b— quark is simulated 



by using nonrelativistic QCD and the Becirevic Kaidalov (BK) [38] parameteriza- 
tion is adopted for the q 2 dependence. Recent results are also available on a fine 
lattice (lattice spacing a ~ 0.04 fm) in the quenched approximations by the QCDSF 
collaboration [39] . 

Combining lattice QCD for the B — > irlvi form factor 33 , 34 with measurements 
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gives | Kb I = (3.43 ± 0.33) x 10~ 3 , while latest Babar data [IT] for 
3.5) x 10 -3 . The results are compatible with the value of \V u i,\ 
luI + v 41 , while a value of |K,J is not extracted from 



from Belle 
\V u b \ range from (3.3 
determined from the B^ 
the B + — > r)l + v decays because the theoretical partial decay rate is not sufficiently 
precise for these decays. 
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In the sum rule approach, the B — > tt matrix element is obtained from the cor- 
relation function of quark currents, such that, at large spacelike external momenta, 
the operator-product expansion (OPE) near the light-cone is applicable. Within 
OPE, the correlation function is factorized in a series of hard-scattering amplitudes 
convoluted with the pion light-cone distribution amplitudes of growing twist. Sev- 
eral recent calculations of the semileptonic form factor have become available, based 



on the light cone QCD sum rules 42-48 . Direct calculations, without extrapola- 
tions, hold in the kinematic region of large recoil, with an upper limit for q 2 varying 
between 6 and 16 GeV 2 . The QCD sum rules provide an approximation for the 
product Jb f+(q 2 ), and therefore the decay constant fs represents a necessary in- 



put for the extraction of the form factor. In Ref. 47 , and in several other papers, 



fs is calculated for consistency within the QCD sum rule approach, that is replac- 



ing fs by its two-point QCD (SVZ) sum rule 49 50 . In this way one expects that 
radiative corrections affect in the same way fs and fs f+(q 2 ), and cancel in the 
ratio (/b/+)//b, together with some theoretical uncertainty from input parame- 
ters in common, e.g., the b mass. Using the form factor results in Ref. [47| and the 
the estimate value is \ V ub \ = (3.46 ± 0.06 ± 0.08±j|]j|) x 10~ 3 
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latest Babar data 

(Ref. [4l]), where the three uncertainties are statistical, systematic and theoretical, 
respectively. 

Last year, Babar performed a simultaneous fit to the data over the full q 2 range 
and the FNAL/MILC lattice QCD results, publishing the following average value 



V u b\ = (2.95 ± 0.31) x 10 51 . This year, they have performed a similar fit, with 



updated data, new lattice results, and some different theoretical assumptions 41 



In Table VII of Ref. 41 , for each lattice and QCD sum rule computation, values of 
\V u b\ are presented, together with the corresponding % 2 /ndf, that can be as low as 
X 2 /ndf= 2.7/4 (with probability equal to 60.1%) in the HPQCD case. The average 
estimate [41] , determined from the simultaneous fit to experimental data and the 
lattice theoretical predictions, yields 

\V ub \ = (3.25 ±0.31) x 10~ 3 (16) 

that is about 1 standard deviation higher than the 2011 estimate. This fairly large 
difference has been explained by the fact the determination of \V u b\ from the com- 
bined data-lattice fit is most sensitive to the points at high q 2 , where the changes 
due to the improved hybrid treatment have lead to differences larger than those 
expected on the basis of the variation in the total branching fraction value. 



5. Inclusive semileptonic B— decays kinematics 

Let us consider the B — > X q I v decays, where the final state X q is an hadronic state 
originated by the quark q. In inclusive decays, we sum over all possible final states 
X q , no matter if single-particle or multi-particle states. Since inclusive decays do 
not depend on the details of final state, quark-hadron duality is generally assumed. 
We can factorize long distance dynamics of the meson using an OPE approach. 
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Another advantage, due to the large mass of the 6-quark, is the possibility of using 
the systematic framework provided by the Heavy Quark Effective Theory (HQET) . 
The aftermath of the previous approaches is that inclusive transition rates have the 
form of a heavy quark expansion, that one can schematically write as 



1 1 \ =i(i2N 



<o 5 > <o 6 > / 1 

c 3 < 3 > +c 5 5 1" c 6 q 1" ° — Z 



(17) 

Here c,j are short distance coefhcients, calculable in perturbation theory as a series 
in the strong coupling a s , and Od denote local operators of (scale) dimension d: 



<0j>s <^> 

2m B 



where mg is the B— meson mass, included in the definition for the normalization 
and dimensional counting. It is basically the same cast of operators, albeit with dif- 
ferent weights, that appears in semileptonic, radiative and nonleptonic rates as well 
as distributions. The short distance coefficients Cd contain the masses of the final 
state quarks (from phase space, etc.), that require definition in a chosen scheme. 
The hadronic expectation values of the operators Od encode the nonperturbative 
corrections. While we can identify these operators and their dimensions, which then 
determine the power of l/m&, in general we cannot compute their hadronic ex- 
pectation values from first principles, and we have to rely on a number of HQET 
parameters, which increases with powers of l/m^. The hadronic expectation value of 
the leading operator O3 = bb incorporates the parton model result which dominates 



asymptotically, i.e. for mi, — > 00. A remarkable feature of Eq.(17) is the absence of 
a contribution of order l/m&, due to the absence of an independent gauge invariant 
operator of dimension four. The fact that nonperturbative, bound state effects in 
inclusive decays are strongly suppressed (at least two powers of the heavy quark 
mass) explains a posteriori the success of parton model in describing such processes. 

The phase space region includes a region of singularity, also called endpoint 
region, corresponding to a kinematic region near the limits of both the lepton energy 
Ei and q 2 phase space, where the rate is dominated by the production of low mass 
final hadronic states. The region of singularity is a cut, since the hadronic final mass 
can vary, in contrast with exclusive decays. Near the cut, and especially near the 
endpoints of the cut, the use of the OPE cannot be rigorously justified because there 
will be propagators that have denominators close to zero. Corrections can be large 
and need to be resummed at all orders. A resummation formalism analogous to the 
one used to factorize Sudakov threshold effects for parton distribution functions in 
usual hard processes, such as deep inelastic scattering or Drell-Yan, can be applied 



to the case of inclusive heavy meson decays (see e.g. 52 55]). 



6. The semi-leptonic B — >■ X c lv 

Schematically, the decay rate takes the form 



?(B ^ XJv) = ^\V cb 



A , , nf 1 \ (19) 



f(p) + k(p)^+g(p)^+0 [ , 



b '"6 

2 /«~.2 



where p — mj./m^ and the coefficients /, k and g, calculable in perturbation theory, 
are expressed as a series in a s . This expansion is valid only for sufficiently inclu- 
sive measurements and away from perturbative singularities, therefore the relevant 
quantities to be measured are global shape parameters (the first few moments of 
various kinematic distributions) and the total rate. While the general structure of 
the expansion is the same for all the above mentioned observables, the perturbative 
coefficients are in general different. 

The leading term is the parton model, which is known completely to order 
a s and a 2 s , for the width and moments of the lepton energy and hadronic mass 
distributions 56-59 . The terms of order a™ +1 /?Q , where /3q is the first coefficient of 



the QCD /3 function, have been included by the usual BLM procedure 60 61 



For the total rate, the kinetic corrections have the same coefficient as the leading 
order, k(p) = —f(p). For other observables, such as partial rates and moments, the 
kinetic corrections can be obtained from the leading-power differential rate, but the 
relations are more complicated and only corrections of order 0(a s p^/m 2 ) have been 



evaluated 62 



For the total rate, the Eq. (19 1 is known up to order 1/mjj, where the terms 
of order 1/m^ with n > 2 have been computed only at tree level 63 64 . Up to 
the higher order available, there are not unnaturally large coefficients, which seem 
to confirm the duality assumption. At order l/mf[, contributions with an infrared 
sensitivity to the charm mass m c start to appear. At higher orders, terms such as 
1/m^m^ and a s (m c )l/m^m c are comparable in size to the contributions of order 
l/m|. The total 0(l/rriQ 5 ) correction to the width is about +1.3 %, in the approx- 
imation that the many HQET non-perturbative parameters are estimated in the 
ground state saturation approximation 64 . On the same premises, the estimated 
effect on the \V c b\ determination is a 0.4% increase. 

As it is well known, any quantitative statement about the value of a quark 
mass must make careful reference to the particular theoretical framework that is 
used to define it. This scheme dependence also affects the determination of the 
HQET non-perturbative parameters. Some traditional schemes for masses are not 
as advantageous in the present context. The pole scheme is most convenient from 
the point of view of computation, but plagued by large misbehaved higher-order 
corrections. The minimal subtraction scheme (MS) sets the scale of order of the b 
quark mass, which is considered unnaturally high, due to the presence of typical 
scales significantly below, down to the order of 1 GeV. Alternative schemes are 
the so called low subtracted mass schemes, where non perturbative contribution to 
the heavy quark pole mass can be subtracted by making contact to some physical 
observable. In recent literature, the latter schemes are commonly used, in particular 
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the kinetic 65|66 and IS scheme 67||68 . Care must be taken in converting from one 
mass scheme to another due to the presence of truncated perturbative expression. 

As mentioned in the previous section, a common feature to many processes in 
QCD is the presence in the perturbative expansion of large double (Sudakov-like) 
logarithms at the threshold, in the OPE singularity region. Resummation of large 
infrared logarithms is essential in order to predict accurate cross sections in many 
phenomenologically relevant processes. For b — ¥ c decays, corrections are expected 
not as singular as in the b — > u case, because they are cut-off by the charm mass. 
Nevertheless, their size needs to be estimated and large corrections to be resummed 
(see e.g. 



69 70 



In order to determine |V c b|, a global fit may be performed to the width and all 
available measurements of moments in B — > X c lvi. A global fit has been recently 
accomplished in both the kinetic and the IS scheme [6]. Each scheme has its own 
non-perturbative parameters that have been estimated together with the charm 
and bottom masses. The fit constrains only a linear combination of rrib and m c , 
not enough to bound the b mass precisely, which reflects on the precision of the 
determination of \V c b\- This limitation has been overcome by including the photon 
energy moments in B — > X s j into the fit, or by applying a precise constraint on 
the c-quark mass. With the former constraint, in the IS scheme, the result of the 
fit gives 



and a very close result 



\V cb \ = (41.96 ± 0.45) x 10~ 3 (20) 



|V cb | = (41.88 ± 0.73) x 10~ 3 (21) 



is obtained also in the kinetic scheme, using the c-quark mass constraint. 

The averages are in agreement with the values given in Sect. [2] extracted from 
exclusive decays. The most precise measurements are from inclusive, that are be- 
low 2%. Still, the determination of \V c b\ from B — > D*W has reached the relative 
precision of about 2%. 

We can also compare with the global fit of the CKM matrix elements within the 
Standard Model, as calculated by the CKMfitter and UTfit groups. The CKMfitter 
group uses a standard x 2 -like statistical frequentist hypothesis, in addition to the 
RFit scheme for the treatment of theoretical systematics. A recent determination 
of | F ch | [71] gives 

\Vcb\ = (40.69 ± 0.99) x 10~ 3 (22) 

with pull, that is the difference between measurement and predictions, equal to 0.2. 
The UTfit uses the Bayesian statistics to extract the observables, and Gaussian 
parton distribution functions to represent statistical and systematic uncertainties. 



Their latest estimate gives 72 



| Ft* | = (42.3 ±0.9) x 10" 3 (23) 
with pull less than 1. All estimates agree within the errors. 
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6.1. Inclusive V u b determination 

The smallest element in the CKM mixing matrix \V u b\ plays a crucial role in the 
study of the unitarity constraints and of the related fundamental questions. In 
principle, the method of extraction of \V u b\ from inclusive B — > X u Wi decays follows 
in the footsteps of the \V c b\ determination. However, experimentally, the copious 
background from the B — > X c lDi process, which has a rate about 50 times higher, 
does not makes feasible a measurement over the full phase space. To overcome this 
background, inclusive B — > X u lv\ measurements utilize restricted regions of phase 
space in which the B — > X c lD[ process is highly suppressed by kinematics. This 
requires knowledge of the fraction of the total B — > X U W[ that lies within the utilized 
section of phase space, which complicates the theoretical issues considerably. The 
region where the background is forbidden overlaps with the end-point, or singularity, 
region, where the energy and mass of the hadron final state are of order Ex — ™& 
and m 2 x ~ ^qcd fty <C m%, respectively. This kinematic region has sufficient phase 
space for many different resonances to be produced in the final state, so an inclusive 
description of the decays is still appropriate. However, here, the differential rate is 
very sensitive to the details of the wave function of the b quark in the B— meson. 
The nonperturbative effects related to a small vibration of the heavy quark in the 
B— meson (the so-called Fermi motion) are enhanced. The parton level differential 
distribution at the end-point region has its own problems, as well, related to the 
presence of soft and collinear singularities. Large logarithms appear, which spoil the 
perturbative expansion and need to be resummed at all orders. A way to remedy 
the previous difficulties is to introduce a non-perturbative form factor, the shape 



function (see, e.g., Refs. 73-79]). It can be interpreted, from a physical point of 
view, as the distribution of the effective mass of the heavy quark at disintegration 
time. At leading order in Agc^/m;,, the shape function can be extracted from 
a reference process, such as the radiative B — > X s 7, and used to predict other 
inclusive B-meson decays. At higher orders, the shape function is no more universal, 
its functional form is unknown and one has to resort to model functions to interpret 
the measured differential distributions. 

The discrepancy between the values of \V u b\ extracted from inclusive and exclu- 
sive decays is a long standing problem, and a lot of effort has already been devoted 
to find possible solutions. Two main routes to progress in the extraction of \V U b\ 
can be identified. The first one is to enlarge the experimental range, so as to reduce, 
on the whole, the weight of the endpoint region. Latest results by Belle 80 access 
~ 90% of the B — > X u lDi phase space, claiming an overall uncertainty of 7% on 
\V u b\- A similar portion of the phase space is covered also by the most recent Babar 



analysis 81 . The second route is to enlarge our theoretical prospective, by compar- 
ing results obtained in several available theoretical schemes. All of them are tailored 
to analyze data in the threshold region, but differ significantly in their treatment 
of perturbative corrections and the parameterization of non-perturbative effects. 



The latest experimental determinations of \ V u b\ come from Babar 81 and HFAG 
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[6] collaborations. They both extract \V u b\ from the partial branching fractions 
relying on at least four different QCD calculations of the partial decay rate: BLNP 
by Bosch, Lange, Neubert, and Paz 82 -84 ; DGE, the dressed gluon exponentiation, 
by Andersen and Gardi 85 ; ADFR by Aglietti, Di Lodovico, Ferrara, and Ricciardi 
86 -88 ; and GGOU by Gambino, Giordano, Ossola and Uraltsev 89 . A comparison 



of the above approaches generally leads to roughly consistent results when the same 
inputs are used and the theoretical errors are taken into account. In Fig[l]we plot 
the values for \V u b\ (without the errors) and each point correspond to a different 
experiment (sensible to a different experimental cut) used as an input in the HFAG 
average [6]. The listed experiments and cuts are respectively, from left to right: 
CLEO, where the cut imposed on the leptonic energy has been 2.1 < E e < 2.6 (Ref. 
[90| ) ; Belle, where the analysis has been performed with limits on both the invariant 
mass (mx < 1.7 GeV) and the lepton pair squared momentum (q 2 > 8GeV 2 ) (Ref. 
[91] ) or with 1.9 < E e < 2.6 (Ref. [92]); Babar, with a constraint 2.0 < E e < 2.6 
(Ref. [93]); Belle, with a cut on the lepton energy at 1 GeV (Ref. [80]); the remaining 
points refer to Babar, where various cuts have been imposed, as can be seen in 
Ref. [81 . 



|v ub | (10- 3 ) 




Fig. 1 . Comparison of different | V u i | values (without the errors) obtained in different theoretical 
approaches. Blu, red, green and violet lines refer, respectively, to GGOU, BLNP, DGE, ADFR. 
Each point correspond to a different experiment (listed in the text). 



Other proposed theoretical approaches 94 96 have not been included in FigJT] 
since they do not provide the same extensive list of \V u b\ determination for each 
experiment as the previous approaches. The latest method advanced [97] aims at 
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providing a global fit to the available data in inclusive B — > X s ~f and B — !• X u I vi 
decays. No estimate of \V u b\ within this approach is available at the moment. 

Notwithstanding the proliferation of theoretical methods and approaches, the 
values of \V u b\ extracted from inclusive decays maintain about two a above the 
values given by exclusive determinations. Also indirect fits prefer a lower value of 
\V u b\. Very recent CKMfitter results give 71 



V ub = (3.42±g;?) x 10" 



(24) 



For Utfit 72 



V ub = (3.62 ±0.14) x 10" 



(25) 



In Fig. [2] we plot a summary table of the inclusive averages from Ref. [6] and 
compare with indirect results, Eq. (24 1 and Eq. (251, and Babar latest exclusive fit, 



Eq. (16). The closest value to the exclusive determination comes from the ADFR 



approach and reads 



V ub = (4.03 ±0.13±g;^) x 10 



(26) 



where the errors quoted correspond to experimental and theoretical uncertainties, 
respectively. 



4,80 
4,40 - 
4,00 - 
3,60 


|V ub | averages [10 3 ] 








1 1 ' 




3,20 












2,80 


Excl Utfit CKMfitter BUMP DGE GGOU ADFR 





Fig. 2. Comparison of \V u b\ values extracted from exclusive, inclusive decays and indirect fits. 
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6.2. B° decays 

B s decays are attracting a lot of attention, due to the avalanche of recent data 
and to the expectation of new ones. Apart from standard form factors and branch- 
ing ratio computations, the inclusive B s decays can be usefully compared with the 
inclusive -B-decays to check quark-hadron duality and evaluate heavy quark expan- 
sion parameters, while exclusive decays represent an interesting probe to analyze 
the (expected sizable) SU(3) breaking. 

There arc well known differences between the B® and the B° system. The mixing 
parameter x s = Am s /T s is about 30 times larger than Xd, and the mass and width 
difference are sizable. Another important difference is that the CP violating mixing 
phase probes the angle f3 s in the unitarity triangle, which is about two order of 
magnitudes smaller than ft in the Standard Model, and hence negligibly small. Any 
large variation due to new physics can produce observable effects, and that alone 
would be enough to motivate the study of CP violation in the B° s system. For a 
very recent brief review on CP violation in the B® system see, e.g., Ref. 



OS 



In 2010 a discrepancy with the standard model has been reported in the mea- 
surement of the (like-sign) dimuon charge asymmetry A b sl of semi-leptonic decays 
of b hadrons. The anomalous A^ , found by the experiment D0 with 6fb _1 of data, 



deviated 3.2a from the SM 



00 



The 2011 D0 update at 9 Ad shows again a 



deviation, at 3.9 a [100], from the Standard Model value 101 

A b sl = [-0.787 ± 0.172 (stat) ± 0.093 (syst) ] % 

^(SM) = [-0.028t°;°° 5 6 ]% (27) 

A b s i is defined as the difference in the number of events with a pair of positive muons 
minus the number with a pair of negative muons divided by the sum. 

A related observable is the semileptonic charge asymmetry a s i, defined as 

= r(5g(f)->/)-r(i?°(t)->/) = i - \g/ P \ 4 
sl r(B,°(t)->/) + r(i?o(t)->/) i + k/p| 4 1 ' 

testing the "wrong" final state, accessible only through mixing. The asymmetry a s i 
measures CP violation in mixing and it is independent from time and from the final 
state (to within a sign), as it can always be ascribed to a property of the decaying 
states. At lowest order in |ri2/Mi2|, we have 

= 1 — a a = Im ( — — J = — — - tan0 s (29) 



M12 I Am s 



where 4> s = axg(-M 12 /T 12 ), Am s = m H - m L = 2\M 12 \ and Ar s = T L - T H = 
2|P 12 | cos0 s . Notice that the symbol <fr s is overloaded, since in literature it is used 
also for the CP violating phases defined in a slight different way. Whatever the 
definition, the phase can be related to f3 s , that in the Standard Model is j3 a = 
arg [— V^ b V ts /V* b V cs ], since the dispersive term of the weak hamiltonian M 12 is 
mainly driven by box diagrams involving virtual top quarks and the absorptive 
term r 12 is dominated by on-shell charmed intermediate states. An addition of a 



15 



non-standard phase, e. g. /3 S (SM) —> (3 S (SM) + /3 S , it is often used to parameterize 

effects of new physics or non-leading hadronic contributions. 

Since it arises from the meson mixing, if there is not a separation of the asym- 
b 

si 

A b sl = C d a d sl + C s a° sl (30) 



metry due to B° and B° s , A b sl can be written as 



where the coefficients depend on mean mixing probability and the production rates 
of B° and B® mesons. Here is the semileptonic charge asymmetry in the B° 
system, which has been measured since 2001 at e + e~ machines. The actual averaged 
value from CLEO, Babar, Belle, Opal and Aleph collaborations is aj? ; = (— 0.10± 
0.37)% ©. 

The value of a s sll semileptonic charge asymmetry in the B® system, extracted 



from the D0 measurement 100 reads 

c& = (-1.81 ± 1.06)% (31) 

The semileptonic charge asymmetry a s sl has also been directly measured by the 
experiment D0 via the decay B® — )• D~/j, + X, using data corres pond ing to 5fb 
of integrated luminosity a s s[ = (— 0.17 ± 0.91(stat)^ ' 23 (syst))% 



-l 
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A new and 

improved measurement of a s sl using the full Tevatron data sample with an integrated 



luminosity of 10.4 fb 1 gives 
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ati = (-1.08 ± 0.72 (stat) ± 0.17 (syst) )% (32) 

The value of the Standard Model prediction for a sl = (1.9 ± 0.3) x 10" 5 [lOl] is 
negligible compared with current experimental precision. The extracted value for 
a s sl is in agreement with the direct and the Standard Model determinations. First 
LHC-b results are already available, from measurements of B a —> D ± y^v and 



B® — > D^lJ^v asymmetries, giving 104 



a s sl = (-0.24 ± 0.54 (stat) ± 0.33 (syst) )% (33) 

Recent precise determinations of the CP asymmetry in several non-leptonic decays 
already severely constrain a possible interpretation of the like-sign dimuon charge 
asymmetry in terms of non-standard CP-violating contributions to B s mixing, while 
more precision in traditional and new channels is expected soon (see, e.g., Refs. 



105 -107 ). Greatly improved precision or, even better, independent measurements 
of semileptonic asymmetries are needed to establish evidence of CP violation due 
to new physics in semi-leptonic B® decays. 
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